{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "c6185b5f",
   "metadata": {},
   "source": [
    "# 混叠效应"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "f1777027",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEGCAYAAAB1iW6ZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABP5klEQVR4nO2deZhcRbn/v2/3bNnXSUIWkkDCEkjYxgCCAiKyqQFFRbxyr8tFvD9UvC7gVbl6FUXlKl5FERV3xYVV9gBh3zIJISRk38g+k22WZJZe6vdHn9Nzzqn3TFed0z0z3fN+nmeemamuOrW99dZbb1WfIqUUBEEQhPIn0d8FEARBEIqDKHRBEIQKQRS6IAhChSAKXRAEoUIQhS4IglAhVPVXxuPHj1czZszor+wFQRDKkiVLluxRStVzn/WbQp8xYwYaGxv7K3tBEISyhIi2hH0mLhdBEIQKQRS6IAhChVC2Cv0vr7yJ9/3seV/YnvYuNHz7caza2eoLv+bPS/G/j63xhT2yYhfe8b9PIZPt+absoe40zrjpSbywYY8v7g33rcBX73ndF/byxr1463efwMGudD4sm1U474dP46HXd/ri/mjhWvzHn5b4wtbubkPDtxeiqa3TF/6B217AH1/yr6h++/wmXH77i76wnS0dOOVbC7Ghud0X/u+/b8RPnljnC7tv2Xac/6Nn4P1WcGtnCqd95wks2bLfF/e6fyzH//zzDV/Ys+ua8bbvP4nOVCYflspkcfYPFuGJVbt9cW96eDU+/9dlvrAV21vwlhsfx/6D3b7w9/zkOfy9casv7BdPb8C//eYVX9jWfYdwyrcWYuu+Q77wK+94Bbc/s8EX9rfGrVjw0+d8YfsOduMtNz6OFdtbfOHX3vkqvvfIal/Ywjd245ybn0Iqk82HdaYyOPN7T+K5dX65+OY/V+L6u5b7wpZs2YfTvvME2jpT+TClFM7/0TO4/7Udvrj/98Q6XPV7v9txfVM7Gr69ELta/HLxoV+8iN+9sNkX9oeXtuCDt/nloqm1Ew3fXoh1u9t84Z/+4xLc8vhaX9iDy3fivB8+jaxnDLR3pXH6d5/A4s37fHG/cvfr+Mb9K31hL2zYgzNuehId3T1ykc5k8Y6bn8KjK3f54v7g0dX47F9e9YWt2tmKhm8/jr3tXb7wS259Hne+8qYv7FfPbsRHf/2yL2z7gdwY2LznoC/8479djJ8/5ZeLu5Zsw8X/96wvrOVQCvNvfByvbT3gC//C317Ddx5a5QtbtLoJZ/1gEbrTPXLRlc7g7d9fhKfWNPnifvuBN/DFv7/mC1u29QDm3/g4WjpSKCVlq9C/cvfrWPrmAV/YotVN2NPehV8+u9EX/sDynfjJk+t9YV/6x2vY2HwQ7Z09Cnnt7nZsP9CBmx72D/Lfv7gFf3rZL2Dfe2Q1drR0+iaPjlQG65ra8YW/+Tvzx0+sw0Ov+wX8189uwp72bjyxyi8Mizfvx9fuXeEL+8Y/38BLG/0D7MHlO7H3YDf+9JK/XAvf2I3/XegfuJ+7cxnW7G6D9y0PS7fsx67WTm2Q/7VxK+54fpMv7Jv/fANb93XgTY9CbWrrwua9h7Sy3vb0Btzz6nZf2K2L1qO5rQsvbNjrC399ewu+9A+/Qvzuw6vx1JpmX9jfG7di78Fu3LV0my/8mbXN+M5D/r768j+W47VtfsX93Po9aG7rws+f9g/ye5ft0Ab+V+5+HZv2HPRNPpv2HMS2/R341gP+ie43z2/GnYv9E9LNj67FrtZOLPeUIZNVWLO7Ddfe6VdoP1y4Fo+94Z8Qf//iZuxp78YjK/xGwcub9uG/Awr16/euwCsBxfvoyl3Y096N3waU/8MrduGWx/0T/ef/tgzrmtrR7Zm8Xt/Wgp0tnfjBo34D6C+vvKk988YHV2H7gQ6sb+oxKg50pLBxz0F85W6/AXTrog3ahPaLpzdgT3sXnl7r7+9lWw/g+kD6bz+4Cs8GJtR7X92OvQe78deAUfDk6iZtov7C31/Dyh1+Q++lTXvR1Nal6Ya7lm7D7c/4dcjX7l2BLXsPYXdrz0S7fX9uTAQnul89twn/WOKX1VseX4umti4s2eLvr2JTtgqdw9VXBLKJ3BPkaDyD1D3Jqfew8PTmebHplXleLlmPRu8pa+EHcO1i1VYRyupL7/w26lc2f5u2dgvrTe8ERezXjJu/UVu7eUWsK9z0FpF9QfH6VenNVzD7yHJh1a9c+pj5w01v0a+RS2tGRSl0WA08PW5evi06yCtOfavk7IU56x3AMQeejeLhyprNMtokLH3MtnKJOvAUp/kKpPdWNpvVgsLTu20VWcm5eZn3i+LkwmIMcc80Sh9TyUWRC8VU1my86vWKMnmVWJ9XlkK3sy50YY5rXeQHrpGSc/PqO2H2W+gWiiOmNHJlzSgLJVms1YxRXk7cyJOXm1dP3AyjDELTW5SVT2+fF7tyM5kQuLwsZIVvaxu5cHMyb60sN96tDECmskZtHU+GTakshR6hg9gZ28LEZpfW5skj97DNck/LE8VUcubpvbllSmyhq4iTV28KxSh7bvLKunJhruSimug2YuXG9Sn0CG3tX82Yp+fautRyEX3ycuP2HlaIqK40UypLoTu/bTqInbFt8iJGSUX0tUZyQxinCBu4EX3oUZbW3tVMFEvMYjDE7deoSq7Xfu0TC91Jb9GvvrYqUr9GHkMlXrnF7Vd/mJs+2oRQCspeoauYHRR9o9CJ6wmzEWYwAy+KMNtICG+dmOTlxPWU1c7/yaxmIkxeNkSfvPQ8oyk5xuVSODm4VWIUN4RNXHYMWfjg2dWIhWD5Jz+DdG5ym+VAMA2i9atiwvpib8SUoil0IppGRIuIaBURrSSizxXr2b3hVQxRGs1rFUdbQullsbEkoyq5nvTRrNZ8egth9j8rnpKzGrgR+tUvF07+5sn96a0mLyeub2/Fvvw+H3wUubDyocdL7y1rzwZwxAkhkg/dnLj9mo2pb0p9yqWY73JJA/iCUmopEY0AsISIFiql3iiUMA68rEdzueRTR1yC2VgnXLGt3BBRrFafMNs/wDsYbMrqwu032KWPZomxBQhNp6/c8skj9muPhW4ul1z6MLJZhUSCnPTx+jWCWPmfFbNfSzV5ubBya9Gv7MovYr+WgqJZ6EqpnUqppc7fbQBWAZhSrOeHwTWwCXwH2StUn3VilT58uVoojZe4G0JmpxGUln/G6iiejtXA5VYIBdJH7td8ml6zD0/PbZbH3OgrtJphJx+ryS9eW2V9cmG/cvNis48URUtmI/ernj6KAVY2LhcvRDQDwEkAXg6EX0VEjUTU2NzczKa1xX9yw95/5+2Unr0rC1+rRyzSWX0wF8jeUqF7/7bxy+aIemTLVSje9DbuJc49YzMYuPQFrdaIlhQvF8Xp16ibqukCGj1qv7Lp4aY379eobRWlX/3p4eTl3dsxN4qifOEs8uQV8+itKUVX6EQ0HMBdAK5VSvm+a6uUul0p1aCUaqivZ1/naw1ndSYsGjjyJkk+f09ZnH8SEb/AUkihF9qhtxFmbqOyEJGVnJPM2y5RrNaExeTHndwwkws9fdZi8mJdJlHkIqKFHvfkh9Veu2vURJzouX61OunF9KuVXDByGZpXvq6cQreQi3Ky0ImoGjll/iel1N3FfHYYcX1arDBH9LVafYEEel6FrU4mvc/l03uekb9UwdXV6ogmtLhWCj2f3Ny9Ffn0Uy91jTzRR7HwmS8mhaaJbaEzE0JEv3KUfoVVXXvvV6uVW1wDLkK/lppinnIhAL8GsEop9cNiPbcQfiUVr4NsNq+CeQJeYTaZUPQutrLEOGG2sfCd3zZfdlFsWxWGW9rG9T/aWehOequ6MhN1wdTIF9ZraebffWXjV47br74TRREmhEIFhXfy08sS9Rx7pH61OD3FG4CF4Sd6N3+TBwR+l4hiWuhnAPgogHcQ0TLn56IiPp9FsUqK/5xLF3dDx780h3l6t6wW59BZ95Dn80JKMvKXKpjVjI0bIv8cz982G8DciaK4k194vm56PS+br/6z+w2Fsy9ev3o+j7TfYFXXaGOIc/tF6Vcv0VYzFnX1lC+KDz3CIR4rinZsUSn1HErv89fwvvqzZ6MypxCISPOF5o93OWHe9167wqJUT3ovmaxCMp8+FzfleT9yxpPemxeXvsf/qHrK6j0WyKRPZbMYgqSvrllPWTOF0md0JWlVV09buXkppQrW1S2rt62CkwNX1poqf3pvWTMF03vkIkJdu9l+VYXTu3Xl5CqkrLxcxOvXrKeshdNzY8Cirky/mowBN30mq4z7NZ1RqE7666os6urtV+9qLFa/msgVk74UlP03Reff+ET+4gO3g37/4pb8O4q9HXzEfz2U/9tt4AtueRbLnBfcu2394sa9uOoP/gspAOBIJv0Vv3oZjzkv83fz2tPehQW3Ps+mdwXKFcbP3bkMf3AutPCWdf53HtfSz/vGY/kLMdyYP120Ht933l3tHQzeurqcc/NTWL2r1Vf+R1bu0i6kCKvr+3/+Ip5d1+zLa0PzQXz4ly/1nt4p7VV/WJK/0MIr12ff/JSW/qivPaxdBnDjQ6vwU+fd1d62Oubrj2jpT//uk/mLD9zy/7VxK/7rntz7272TJ1fXd//kOTQ67xp34y598wA+/tvFRnW98o5X8heduGVt60zjosAlC25695IIt1Rfvms57nh+s1bXk761UEt/8rcWYmdLh6/8v3x2E258MHdJQ7aAXLzzh8/kL/9w4y5a04xrAhdShNX1Q7e/hEWrm5yy5j7btr8Dl932Aps+bww4YZ/+09L8e+W9dT3ze09q6Y+94RHsc99V70S9+bG1+JHzrvdsyHh3edv3F+Xf3+6uAu9dtgNfDryXX6urE3fBrc/jpY17fWVduaMVH/31K72nd35/7LeLcd+y7VrcYlH2Ch0A1jf7OwgAfveiriTDWLxpnxZ3YeDigSDeQfK4c2uPN+z1wO04Lp3p3MD1Lt3cCyG8CnlPu/92H5et+9yB2xP3j05dTU4IuLezeMt677IdIbFzeJ/6tHP5hLetXt7Ev7Sfe6PlA8tzSs57FO/NwE1ELnucm2y8deUGfneI1bN6V5tW/r84N+GkQ9rKW9bgwAVyiq43vP3q3trj7Ve3TEHanZuvvP1yl3NJgjcs7MabTfnJqyfu7y3GwNI39wPwt8uDy3eGRXfy6vn7SUehe/s1eAGNi3cV4OIqOW9ZdwRubXJxb3Py1urPL2/Ryh/Gyh0tWvq/u21tIBfPr89dtOHt1+fW7wkm8eHtw4cDl90Uk4pQ6L35eUzOtborpSjffszln3uAycDhyuqGmShk1l/nhJnlH17W8P0GPX+TvLyukmB6k6/+u1XlcrLpV87ZatLX7hI6yrcfgfj9atPW7IZvhLaO+o4dm371ukZ78rcYQ653SvmfANj1KydYRnLl/Lb6EpQv/0jJjKgMhU5+P5UXGwGNduTKK8w26Xvi2igOTkflB6OJ4ulFSYQXXy+rycDJ+56ZA2pmClnvV5u27nVCCLPE2A3YqHJhP9GzdbWYvKLLRXi/mnx71CYvr+85kH30ulpNfuF5hcuFXgAzoya8rqWgIhS6S9zBaGWhR5zd3eLEnXyinPzwl8NCmNmymudRirrGfVtjWF/ZxOXTM2W1Sc+E2eyj2UxeHFxeYcn5o7cWCp3NP25bWdTVYuXG1tVmvBuXKh7FfDlXn/Ff97yuhb20ca92YXA2q/D27y/S4t4cuAAXANbsasPn7lymhQdvkAeAXz+3CRsDN43vbOlgN0Y+FrjBXimFe17dpl0O3NqZwsX/p+cVvHAaAJ5e24xbF/kvN05lsjj9pie0uDc+qL8bbfm2A7juLr0Nz7/lGS3sZ0+tx96Dfn/+lr0H8e+B2+oB4IrA5qhSwF8Xv6n52Pcd7Mb7f65vmF3z56Va2GMrd+GPgYuwO1MZdiP1hvtWaGGNm/XLlZVSOPsHevofLVyr+WDXN7Xhmj+/qsUNll8phd+/uEXzke9u7cQVv/TfVg8An/ydv/0UgPtf24F/Bi5Sbu9Ks/1yHbOJ9/z6Pfhh4ILwdCaLM5jNxeBF6EDOtxy8rR4ALvqxvpF7+zMbNB/31n2H8DFm0/ijv/bXP6tyF38/E7gc+sChbvYwwefu1Nv/ydW78evnNvnCutIZvI0Z78FLnAHg1Tf346v36PJy7v8+rYX95Il1OOhsWrtsbG7H1X/U5fWDv3jR939WKfzlpTfze1elpiwt9D+/7B/gBOCHj63V4nWmM2hzNpu8/HSR/5ZvIsLvXtzM5hW8QR6Advs7EfDICn6jI7iJllXAfwaUNAFYvpXfRA3edE9EmpImIuw72M1al7981i/0BOAXT2/UI6JnY83L9x/xT34E4O6l/C79Cxv2+v7PKqVNHESEFzbwG0gPBDbhCMAN9/kHIxGwZS+/iepuAnrz+vET67R46azqOSnhIRiXCPhTQNZclmzZ7/s/q6BNHATgiVVNbHp3I70nvcKX/xGUC8Kqnf6b6l2CN90TdCVNlJsQvEf1XG572m8QEIBfB2TFZc1ufSM3aDwREf65nN9cf3adv7+zSuFLgQmJCFi82d+mLvcFNu0JhG/+MzAGAOw8wG+i/vaFzVpetwZ0AJAzAHcyG7H/G5gkCcDfGrdp8QDglYDxklUKX7vXP3GU8hW6ZanQgxDxy3Lv+dxe04NfanPw/kQyTp87HxxIHVJ+Dq6sROblJ6J4LgwyXz5yj7Rpa65dCHz52fQheXEKjk9v169cAUzLmjvLHEges18J1Mu+CJMXW654ri02HtP8ubY27FemXcLKz6YP6dew01Js/hb9yhSgZFSEQgd05UFkPnABc39g2AAx9b+HxbMZeNwzCr2RL58evK/dTslFr2tY+TmU4vs17l5BXLmw8aea9is30ff23CBh7WLTr1xexkYRmSv/cLkwSs4+I6z8NmUwVeiA+eQV9eRcVCpCoYcpGdNvZYVZQlwHR734wIU35MyVJKBbIgTzgZdL3zfCzFliNoQpOWMLPaRfreSCCeePfDLpYW5iZxW3GjFva4CRC7L7ZqJpW5kcb+2NcCUXbzVgegAirF9TMY0afqLn0peOylDoTAcRbKwTcyXHXyYQU5gtl4vByERkXtewycsmveFpIn7g2rgxgGBlCbmvf5sQVtYuG7mIOdEbu6eyzOQV4h7jyJVVbyvjuoaUlZOLUKOESW+q5MKeaxqPqDT9yqaPbQCKD70g+tKcrKxOTshMhTmXPp7LxsZCj+NGCFvacm3F3pAE8/TxXS6qNP1qld6srUNdcYZWY869xMmQUfL4EzXM3RA2/cp9c5O/ecvc359VinG5lKZfw9PrYay+KO2rWzQqQqETwEqzuZLjrUZuuWljiaWY0cgNcIL5WeOs0m02W2HmSsstN20sKX7gm6fncF88FsRmacx1TFyXi9XKjUnPf6mL71crl0lwNUNk9zIoQyVl06/8GArL3twHz+ZlYdRwGMsFzA9hhMlFqagIhc6NWxslF7Y0NxXm0F1z0+WqxQZu2CmZuCc37NwIZm0VZuFzwzbMZcP1a5epkrIoK5/col+5kxsWVnNYv5q6EdhTMijCGIg5eYXVlUtv+joKpfjJy9yo4fuVa+uwe3xt+pVLXyoqQ6Fzwmx5ysXUOrH5yrhpBwNAdyajhYX5H+Nsatr4SvmBxz+AbasQYeYsIa78GdavbHtKxaysvTzBKH3oNwyZMBur17Rfc3JhlheHzeQVqqQMV8ncGAqbvLi4rFyE5BWGuVFi/Eirfi0VFaHQs4z/kWC33DT1H9rs8JsuwQAglTZbwqmYy033GUH4uuppQ91LhkvrsH7hT1Nw/Wrnhojbr5wP26ZfM8wDbPLn+jXUZRPnKJ7FiSC2X4nfFDWWq5C8wlYI3HFWK7mI2a9x5apUlJ1CT4d0cNwNoTiWFBF/8sLG6uXcCMZWK+w2Ck3rGv5uF8O6hqTn4oa2lbEri4loYbWG9aupJceuRkCskmD7NcQNwcVllWyWKantaibmasR8AzmuXABcz9qd6InRryF7buZyUTrKTqGHzdimViM78EM3Rc1mbDvrgs2etcRMNyrDrBO71YTNwNXDTC2xsLhhbW3ar6HpY66cTPcbbPZG2H5lffD8aiZUrhgfetxz6GxdQx7JnWixWqUa1jV0lcrkZfPCubBVIpue3RQ1X82UirJT6JxrIrchpEszN2OH77Dr2PnALaxO0/ShwmxmnXADLHRpbLwBHDLwLHytNpaYdrY65Mx9mJKL268lkQuLydNmNWNc1pCTVnHqGva9D/PVUEi/hsmFlt48fxXSVry+iGfURH1nelTKTqF3MZuHNoOBU0bhZ3C5jUo9Hyslw1pidoNB8x+GxA0VZsOlcZiFb+UyMU0fsilsOvDCXTZm/RpmIXMWns3Ky7xfmPSwq6u232Ahl7kylMBlYrjZH5Y+fOXFlVXPK3RCMVXIJerXUlF2Cj1cyQWEGbwwcAOECMiwPnDz5aKxrzPU12puyelHtsLcELyAmR9b1OMh5HgYP/DY5CEWPp/etF/DXE5cGbhVHidXoT5s4/0Cu37VUsfs17D0Yb5mtq0M+xWW/cokR7fFPpRu1JjvV3CTH8D3K7dnZztemQeUjLJT6Ka+0twsqs/YXSk9DLBZGvPlijNj26S32QAOW0Kabl6FrWbiWGIEG5eJ+Zn78KW5PvlxG9CdjKyEPddqv6AUq5lQl4uZAcONAYUwH7bZiti2X7X0llav7oqzG8Ose8li8ovbr6Wi7BQ612gtHSmtg7rSWe2l9AB/+XJ7ZxpdzIBu69Qv5N3PvEe7ozuDDmaQcBf6cmEdqQwOMWU9cIhPH5xUutJZHGTe++5eshxMz7VhO5O+uU1Pf6g7jc6Unr61Q0/P1fVQd5qta1hbBfu1M5XFoW49r/2H9H5p6UhprhSlcv0dhGurg11pdkBzcnGgg5OLdCy5yNWVkYuQ9MG6dqWyONhlNgZaO1IhdY0rF0xbMXJ9qDuDjpSZDB04xMlFhh0D3Hht6UixrjROLsLrqrcrV1eu/Oub2tHUxr+7PS7U1+ckXRoaGlRjo37rTSFe39aC9zC3CAmCIJQLV591JK6/8JhIaYloiVKqgfus/Cx0q3eWCIIgDDySJdK85afQrb62LQiCMPBIlugsetEUOhFdQERriGg9EV1frOcGEQtdEIRyJ5EYwAqdiJIAbgVwIYA5AD5MRHOK8ewgNu8sEQRBGIgMdAt9PoD1SqmNSqluAHcCWFCkZ/sQC10QhHJnQFvoAKYA2Or5f5sT5oOIriKiRiJqbG5ujpTRzPHDMGPc0GilFARBGAAkB7hC50qnnYdUSt2ulGpQSjXU19dHyujYw0bi+guPjZRWEARhIDDQXS7bAEzz/D8VwI4iPVujRJObIAhCnzDQXS6LAcwmoplEVAPgcgD3F+nZGolS3uEkCIJQYpIlUmFVxXiIUipNRNcAeBRAEsAdSqmVxXg2R7JUrSEIgtAHDHQfOpRSDymljlJKHamUurFYz+U4c9Z4vP0o3Qc/f+ZYLWzyqDrUVOnVfOexE7SwE6aN1sJGD63GmKHVWvgFx03Swo49bCRb3pnjh2lh5xytl/+Iej0eALxlxhgt7K1HjtPCpowewqY/i2mrU6brzxw/vBbDa/U5/sLj9boeP4Wv6+RRdVrYucfobX3UxOFseu65Z84ar4VND9kYf9tsPe78GbpcTBpZh7pqXS7OmzNRCzth6igtbFhNEuOG1Wjh72LSHzNpBFvWI5n+5vqKkx8AOO0IvV5c2JTRQ1g3JSeDJx8+WgsbP7wWI+t0ueDGwHGTebmYOkaXzXcwcjFrAi8X3Ng8Y5Y+BqaN5ccAJxcNzBiYOLIWQ2uSWvj5x+n9OneKLhc1yQQmjKjVwt95rD/9QHe59CnVyQRuvOR4LfyK+YdrYbMnjsAN79aPxH9zgZ7+vSdM1pT/hBG1+O775mpxv3j+0VrYO46p1wSaCPje++dpcT999iwt7NSZY3HJiZO1cC7/K0+froUdN3kkvsSU68ZL9bq+72TtEBJmjBuKbzPt+l8X6ZvQ58+ZpA2e0UOrcRNT18+eO1sLe+uR49kBfdP79PSfOHOmFnbC1NH41NuP0NMz+X/wLdO0sFkThuOb7z1OC/8GE3bxvMMwIqDQxofIxRfepbf/WUfXs4ry+5fpZf3UWXqdGqaPwYca9Dp8l2mrfzlNl4tjJo1g+/BbTF9fcpIuF1NG1+HGS/W6cu8ieeexEzWlPLQmie8z/fL/ztHHwOlHjGMNiO+9X8//396qy8W8KaPxmXfoz+Xk4gMNU7WwmeOHsXJxw3v0sAvnTsL44X7lPW54DW5iyvr58/xjYKBvivY53AzHWeI1VQn2VZl1IXG1fIjYK6NqufRJfWZPELHWEVvWkBc8cPmH1ZWjrlovF5dXrq56erauVQn91bbg9ze4ctVWJdjX65qmD+vXKgu5IOZwFisXIW1lXFeLfg1raw5Wrpi8uL4CzOUCFnWtYfrVagyFlNUufz19WLsEqU4m2LxM5YJg1q9ioQewEWYOTphrkwn2RclxhCmsg0MHHlNWNn9m8uAGE5ATUi6uBtnVVUtuM3mFKSkm2CY9NyFxbV2VIPawLavkqpLsfZ02ZY3dr8wT4vQVYD555CZqPX3YRB0nfezJK2Ty5HzWoQYY84haTl9UJRAUDAqb/AJ5iYUegGu0akYYkkTs7QOckquuIuZCBF5JsemT+n2dRGHCzJQ/GWadcPnzk4Rpet7qDGlXtq56Xgkyn7yqk+ZKjp2QkjYWvh5GIZdocgqJ61eE1JXrF6t+Zcoa1q+mk1eurfUHVDGai2vrsH4NlQutoLxFGp6euwDDrF+5MQyYj6EwPWvar2HjPdivA35TtK8xFWYi/naSsBm72G4EgtmM7T7T2JILWS7GSU/gXS7mlpSdhR5n8rFamjNtTeAHflVoXvpVeFxetXH7Ne7KLcSoMZ4QLFZeoRY+O4bM8qqN2a/cGAbM3ZYE/vYp036hsFVuIL24XAIYW2LgO8hmxuYEP2zG5grAryYM0yMsf96S4jC1eon3QvB+acZCJmtLTMfYEuOsZtj1q+ndLryhELJyK0G/hqW36dc4E0rY5GXarxTiQ+fHkEW/xpSLsH7lMF2hhBlFQe+BuFwCmM7YYV9CCpuxg8IQuvll6JdOhEwIXHqbjULWVxnbEuMHvo11E99XatavnLsAMC+rzZfTeLmws7CDHRNmydltFPJlDRJWVZv03DPYVW7IGDBezVhMXjYK2XQ1k7CZ6LmVm6FciMslgKkwgvjLijnChJnrIHMlZ3dCgMM0fRg2LhdTahhLLOw0Azt5hfrAmbzCJp8YSg4hkxeHjXuKc9mEbxSau9I4TGUw7Jp5G7mwWc3oSs7GFRfmMjFLHz55xWsrjnCjRhS6NaZLuJyv1Iyw5aYpYctdzpjk4oZZnabpw4SZ3y8wT89RzZ3oIX4pyS6tmeUyELa0D+tXw03RkH41VVJhbW26bObTE3sNmY0rzaZf2cnP0D1mN/kxBwNg3q+5sOL3q6l7y2YMhK0QuCYMtrVcQRfAdMa3mQj5JZTl0lw7+WHusknY+DpD3AimO/xhm1+x/MohZeWtdn25nYvL5GXRrzZtZQp7HDVkNcKmZ1czdu6hOP1qNQYsJhQ+PbMBTearGZv9Apt+tXK5sE/QCdsUNZGLUr2PqmwVuqn/zmrGDekgU2qZI1dhSo6zjog3TiyWxmHWiYV1YyjOYROCaXuFLTnj9qvp2XBbS4zrVxtXnn56ip/ouQ3oMLkw71d+QggrK5felLC2NnUxFKNfuboar2ZA5v3KHVuE2apeXC4BzJdg5lYne64U5kqOc0PkTrmY5R/uRtDj2hzR5ODO7Nue/DB1eXCE+8CZZTgzmsNWE8YnP2A3eXHuJdO2rkoy/RoiF5ySDtHnISc/4vcrn95wDCT1fgkzajiK0a+8UaSH2RwM4KgO2S8wQU65BDBemlvUMOwr3qZw1lWYy4UjLJqpJWblRoibPqStTR8RfvqIiRtmtRo+N7bLJWZbcYM37OQHh9XJjVL0axHSmxo1Nq40fjUTJld91VZmaeUcegBTYYbFpih/FNDcusn5OvX08a0TJow172Bc2FCXjVHqMEvMXEjDLGTjtgo9uaGH8e8nsTyeFkxukZ7r17BTMnx63kI23Wi0Ia7VWpVg3FOGfmUg/BCD8WrG0gfOFMD4AbWcXBjqG3G5BOCEmWsku+VimP/RDC4vgrnVGqaPbZRcrBM9Fm2VG7hcerP8bVwu4en1B9j5lc2oZt7xY+Oy4epqNdGHlJVLHqrkLFwmHOZGDVNXs6S5uBb9yqa38IHHdcWxcmE4BmRTNICN/9GUuOfYOavRZqPQxo1gk54j7hlcVkkhnpIDbPYbzC2x8A1o880v/gFmcDHt/Mp8uNWEYLEa0fM3b21+8jLP32a/wDT/MPhvEJv3axwjWyz0AKYNfxhz4UIY3Ix9+Fj+IgWOobVJTRiPrB9mXNYxQ2s0hVg/otZYd0wYUWsszJxAcZcQhMG5p8IuYuAYWadfGmLTVsH3UAO8TzUM7iKOMLhNzeljhxrPKMPrqrSosyYMN1YIY4bWaP06aki1cfpJo/S2CoObvKZZyAX3tsqwi1s4RjGXyUwbO8R48qpnLpcIg5O1yaPrYrniZowbBhPBkE3RAvzl30/Tws49ZgKuPutIrYMWffFsLe4tHzpRs04apo/B/yw4TuueR659m5b+uguOwZH1w319WT+iFr+8skGz8O/7f2do6f/19Ok4++h6rawPfOZMTZj/9qnTtfQXzZ2EK0+foSmeZ750jhb39o+eooWdMWscrr/wGC3/hZ9/uxb3m+89DhNG1sFb2SPGD8OPP3wSW/4gnzrrCDTMGOsra1WCcOdVp2tK6q5P63V930lTcOlJU7S8nvrS2VrcP33yVC3s7KPr8ZlzZ2v9+sQXztLi/uCyeRhWW+XLa97UUfju++Zq6R/6rC4X/3neUThm0kifDIweWo3ffGy+Fvee/3irFnbFqYfj/OMmsjIYlIs7r9LHwPnHTcTHz5ippX+KGQO3XnGypuROO2Isvv7uOVpbP3qtLhdfu/hYTB3jN4CmjhmCn374ZASV3D+v0eXiE2fOxOlHjNPyuuvqt2p1/cfVuly854TJ+FDDNK2uz35ZHwO/+dhbtLC3zR6PL5x3tJb+8f/U5eI7l87F6KE1vn499rCR+MEH5mnlf/Czel1tDmvYUDEK/XTmSrYFJ01BFXO8jr0Sjrk954LjJ2FoTZU24R4zSb9m6wLnlhVv1HOPmeB0uj8ud53We0+cnFuaesJOmT4GE0fWaUqOu2rvkhOnIJnwL21H1lXhcOaqNu76vvfMm4y66qTWVrMn6tenude0efM6b85E1uo+nrmm6z3zJmvp335UPepH1GoD95Tpel0XnDQFiYTfap42doimTADgDOb6ugUnTmaPqB5Zr19/dq5zdZg36vnHTdKUPADMYa5fu2iufvvOOUdPwNhhNZriOOnwMXpZT3DkwpPZ3CmjcNioIdrK7bQj9DFwyYnOGPBkVledwAxmDJzNXEn37rxc+DmauVbvfOdKOm9Z33nsRIwaWq211VzmWr/3nKCPgTNnjceEkXVaXRuYawUvPWlyTi48D5g4shbTmFX222frdXVvLAuWlbsW751zcvrCG/VdcyZiRF211lbHTdbrKi6XCNgsa+I2cNwlVG9LSpPlpk35ubg2x6jYvCyqz36btJfPtPzZ9Da+T4u6lrBfjfLvpV+MvpHI7R+EtFXYoQJT4h7F662tTdrRpq25olqNobj9Ki4Xe2z2slglYSMg8fbN+MFk8Rx+4y9k4HInhKwGQzyF2pviMKurHha3rfsyL5uhbHMGn80/5NQLm1fcfrXIi82/l7Y26TK+rfmE3NiIaxT15eQX+tySPHWA4Aqo96VXYc3oRvFuCuWtRoPGdzvYv1mYCzNSUk4krx+/R8mZW63ejd3QgduLMCYNnHtuXb0bYDYK2S2i/1ove8VRzfSVTXqTAezGGeqtq5ObiR/U7VdfXcn3y6isPrkwSJdP75S/ynOiIyx9b0rKyEJ2ilhXo7eVkQyzYyi8bFocVy48L36zUrJ5fWG+8hnCyYXJGBALXee/3zNHC7vbs7HkCvEVpx6Osxi/sfcmcreB/+nZxHM75sxZ4/GBU/Qbwm+94mQt/S+vbMD44TW+9FNGD8Hnzp2tpfduzLgC+z/vPQ7HOP5J7yD42sX6re3ezVE3/afOOgKnOj527yD8DnNrO7dheeHxk3DxvMO08B9ffqKW1x8/eSqG1vhPNcw5bCQ+eeZMLf0vr2zwpM+J3c2XzdP8k3XVSVx3gX6b/J89m5tuv37+vKNw0uGjtbj/s0C/of1ez0a0W/4PNEzFO5i9k5s/cEL+b7cNvXLl9utpM8fhw/MP19L//CMeuXAi/+wjJ+dPXLnPrB9Riy+cd5SW/ncfn6+l//rFczDnMMdH7+nXr79bHwPeDUNX8Xz8jJk4Y5buY7/pfT1y4daL29x957ETcMmJk7Xw//vwSfm/XcPpN//2Fowa4t9PObJ+GD599pFa+l//a49cuIr0O5fOxVEThztl6pkQvnKhLhfewxBuW33mHbNxyvQxvvQA8K1LjtfS33+NLheXnDQF73L2ibz88IM9cuGO979+6vR8ud2sTpo2BleePl1Lf9u/9BxGEB86w8fO0BXHyYePwWlH+BVaXXUSX393TiF6leQH3zIt/7fbwEdPGoGPnjbdFzeZIHzlIl2hehWfO3CmjR2KzzrK22thX/tOXaGfc/SE/IrAzX/c8Np8Xt4u/+TbjtDSz585Nj/I3boOranKD3Jv+itO1RXP8VNGYYEzSN2yVicT+CpT1wUnTsn/7ZZ15vhh+A9nkPa4hwhfPP9oLf15ngHiDoYJI+vwxXcd7csfADvw3zprPKaMHuKr68i66vxE5+3XK0+foaU/cdro/KTeYwkm2YnyMs/k7cadPXFEfqLKW60JwvXM5HPh3MO09FNGD8nLgLetPsNM9GcdVY8RtVW+uo4ZVoOvXqzLxSeYybNhxtj8ROe29ZCaJL528Zx8vi6XeyYkN3zO5JH4YMNUX15VyQT+i5GL957Qo+TdvKaPG4ZrzpnlPLPn2dzk5W46e+taP6IW1zvK2ysXnzpLl4vTjxyHI5wNXreth9dWsf3qjmsv86aOzitvV8fWVCXybeXlfSfrcjFrwnB86qzc2PSu6L/MyIV7cMKNUwrKWqGH4e5Se2dBV1jDmtEr5O7pCa8wFWp/7xLKzd/rvwtbcrp5eS2JRH4Q9J5nLr2Tv6+uFumZshZaXrN19bVV7+n9rgqnrQuWtOf0hE2/6jn5B5NVXd08vf1aYAQlCrRVGMXqV29bJ2zaiilrIbeJNy/F9GvBtvbWFcXpV5u2BvS2DqNQv9roi2JSFIVORD8gotVEtJyI7iGi0cV4blTyg4EidnBeyfVQWJi9HaxPCIXy4oW58AN6lJynLPmMDdK7MS2E0T9w3Zz0CSmMZKL3wRBGfuAy/Woy8vP94gkrPPnEm7x8dXV+x+1XGyXnHwP5AhinL2a/WinJvFFl0VZxx7u3XwtUlu1XTi5DGOjn0BcCOF4pNQ/AWgBfKdJzo5G3Thir10LJgRsMIfADtzDcwGOyLwin5MyE2f69GXEt9EJWbxjc5NfTr+Z462fTxqavNvZSqK3C83Kfz/WruZaKvJphhLigkuL6lfT8w2AnhIIlLbCasXiCf6LvPYX3c9MbkrwMaAtdKfWYUirt/PsSAH0HsQ/pdbkX0UKP5oYwty7Y5WpEC99GyfHWRe9p/JOX3tYFLbHIFrqupGyUHNuvgQ0to/QRLbG4/epi06+sUWPSVjFdJjYKudf0Fm0V3ahx8zLv18IuWvO2KialMPw/DuBh7gMiuoqIGomosbm5uQRZ5+A6yEKfAxE6iHNDmOcUNhjMB17CV1dzYUaEyYuYyQtM/mEkCwyG0KIyA9emX3vbG7Gy8L1/20x+PUu/gnlw/RrJlRdRyYEbQwW0BdevUb9wxvm1w2An+oR59tyK2sqAjrAPVapN0SrTiET0OAD9e8zAV5VS9zlxvgogDeBP3DOUUrcDuB0AGhoabPSeFbxfOMKMbbFJkiziwOPKXzB9ZPeS/eTlT+/mZQ5roVu4wngLvXC+vW0A26ymrCb6IvarlXuql36NquQi7TdE/MJZtH0oT1lK3K++9M5vm32oUrlcjBW6UuqdvX1ORP8K4N0AzlWm75stEb1tftn4auMuoay+js4IswmlUHKFLDHuATbyySkpG1dY1K/+97ZRVyoll2ScrTZ5sRZ+TKMkupIrnG8QK7lg91YKw7tcLNK7eUVV6Ey/2uwXFBNjhd4bRHQBgOsAnKWUOlSMZxYDdmlt0Y7+gVsoL0bJRBTm3p4ZpPfJqzBRJi9feud31K+Ic+UvhM+9ZeFrdWHbyia9zX5DVLlglFQ+vUnyXk5u9JeSK0TcfvUrSQuNnk9h3q9seos0A/2Ci58CGAFgIREtI6LbivTcSHD2rd15Ze40g40wG0fteT5ntVrAWZ1W6T1/R1kOWgkz+84OczUT1bphb4OP2a9Wx1mNc+qB24COSjSjRrfwbYiukO3xj1GLVS4T1cbHHWW8D2gLXSk1qxjPKRbc5pNNA0Y5l+pLb/FlGRf/0hzG6XuTJZvTDP50Bhm76aMIc0TrhN9viDBwGaszql/ahiib3b62sjoRpBN38oqbvhCcQrY5s+9rqihjqA9XqSXS5xX6TVHnN+9yMRlMTlwrlexJH2W5GPHkh5UTmkse2+USb/KKMvD4Y4/R+jW/X2Dll45GlI2+Ql/iipaXRVv1pSsu5hfO/EcJ9bDQ9KzbsnC+PendvMzT2LzJ1YaKVOjc5pPVV4mjHN3wpneTW3Qa6z+0yotJb2XJRfMfxhVmu2OL4ZOHUfa9TF52c2dUlw+s84qq5Fy8CtXulE1xVm6xN0UjTrRRJl+bc+i+/GO6wopJRSr0XhWqkZKztzp96SMIU9wvJvkXq/aWYPwjWzFXCBHzivLFJP9Eb5E+toXu5mXer/zX4W3yYsIKJ4/0TU9/+iibotzKLZorLFJbecJ6TgRZpC+R1W1DZSp0bsa2EbCYHRTlSxVxXx3AKjmLJ3DCbELcA6p2Ss6J6w2LpOT0yFZtFXPlZhM37hfOoiq5CCLsT25jYrPZmxeAa5e4bWWzcnPpf3VeoQrdvS9xeF3Pnm+V862Do5g7MoNMdl7TOm5YjfZZ8LZ5brO13okzaWThm+XdG9G9Sma0c/M5dxeilt6pq/dF++7FE0ca3Lbu3sM5xlNXtywmt7VPHJmr6wTmtvVg2wTbDgDGOvm6r8btjdnOu9OrPRcgDHNeMztzvH7vY5Dp43L1GeGRC7eMRzF3ZAY5bJQrF3o9Rgduqx9SndTiuO/Jnzy6sFwc7cipd+U20nnHOHdPbBD33txhNT11dV/VzN2RGWTqmFxdxzJjwP2sNyY4cuHKR29wfT9maI1xXm59vBdjuO/pN5Hhw51x5n2Hu6vQj2buDw4yyXnP/XhmDATvCwi+J77oKKX65eeUU05RxWDVzha1ZlerL6y9M6Uef2OXFvfpNU1q/8EuX9jG5nb1+rYDvrDudEY9/PoOlc1mfeEvbtijdrd0+MK27T+kGjfv9YVls1n14PIdKpPxp1+6ZZ96c+9BX1hzW6d6fn2zVtZHVuxUnam0L2zF9gNq3e42X1hLR7d6cvVuLf2Tq3erlo5uX9i63W1q5fYWX1hHd1o9smKnlv75dc2qua3TF/bm3oNq6ZZ9vrB0JlfXYFs1bt6rtu8/5Avb3dKhXtqwxxeWzWbVw6/vUN3pjC98+dYDamNzuy9s/8Eu9fSaJq2sj7+xS7V3pnxha3a1qtU7/XJxsCulFq7U5eKZtU1qX7tfLjY1t6vXtu73haXSGfUQU9eXN+5VuwJysePAIbV4Ey8X6YBcvPrmfrVlj18u9rR1qufW6XLx6IqdqqPbLxdv7GhR63b769rWmVJPrtLlYtHq3erAIb9crG9q08ZAVyqjHn5dl4sX1u9RTa26XCwJyEUmw4+Bxs371NZ9/rrubu1QL6z3y4VSSj38+g7VlfLLxevbDqgNTf4xcOBgt3qKkYsnVu1SbQG5WLurVb2xQx8DjzJj4Ll1zWpvQC627Dmolr253xcWNgZe2bRX7TjgHwM7D3Solzf65cIWAI0qRK+S6qcvdTY0NKjGxsZ+yVsQBKFcIaIlSqkG7rOKdLkIgiAMRkShC4IgVAj95nIhomYAWyImHw9gTxGLM1CQepUPlVgnoDLrVWl1mq6U0m+9Rz8q9DgQUWOYD6mckXqVD5VYJ6Ay61WJdQpDXC6CIAgVgih0QRCECqFcFfrt/V2AEiH1Kh8qsU5AZdarEuvEUpY+dEEQBEGnXC10QRAEIYAodEEQhAphQCt0IrqAiNYQ0Xoiup75nIjo/5zPlxPRyf1RTlsM6nUMEb1IRF1E9MX+KKMtBnX6iNNHy4noBSI6oT/KaYtBvRY4dVpGRI1EdGZ/lNOGQnXyxHsLEWWI6LK+LF9UDPrqbCJqcfpqGRHd0B/lLClhL3np7x8ASQAbABwBoAbAawDmBOJcBOBh5N5ceRqAl/u73EWq1wQAbwFwI4Av9neZi1SntwIY4/x9YQX11XD07EXNA7C6v8sdt06eeE8CeAjAZf1d7iL11dkAHujvspbyZyBb6PMBrFdKbVRKdQO4E8CCQJwFAH6vcrwEYDQRHdbXBbWkYL2UUk1KqcUAUv1RwAiY1OkFpdR+59+XAEzt4zJGwaRe7crRFgCGIdpd0H2JybgCgM8AuAtAU18WLgam9apoBrJCnwJgq+f/bU6YbZyBRjmWuRC2dfoEciurgY5RvYjoUiJaDeBBAB/vo7JFpWCdiGgKgEsB3NaH5YqLqQyeTkSvEdHDRHRc3xSt7xjICp27ACRo/ZjEGWiUY5kLYVwnIjoHOYV+XUlLVByM6qWUukcpdQyASwB8q9SFiolJnW4BcJ1SKlP64hQNk3otRe49KCcA+AmAe0tdqL5mICv0bQCmef6fCmBHhDgDjXIscyGM6kRE8wD8CsACpdTePipbHKz6Sin1DIAjiWh8qQsWA5M6NQC4k4g2A7gMwM+I6JI+KV10CtZLKdWqlGp3/n4IQPUA7ytrBrJCXwxgNhHNJKIaAJcDuD8Q534AVzqnXU4D0KKU2tnXBbXEpF7lRsE6EdHhAO4G8FGl1Np+KGMUTOo1i5w7+5xTVjUABvJkVbBOSqmZSqkZSqkZAP4B4D+UUvf2eUntMOmrSZ6+mo+c/hvIfWVNVeEo/YNSKk1E1wB4FLkd7DuUUiuJ6Grn89uQ24G/CMB6AIcAfKy/ymuKSb2IaBKARgAjAWSJ6Frkduxb+6vcvWHYVzcAGIectQcAaTXA34BnWK/3I2dUpAB0APiQZ5N0wGFYp7LDsF6XAfg0EaWR66vLB3JfRUG++i8IglAhDGSXiyAIgmCBKHRBEIQKQRS6IAhChdBvm6Ljx49XM2bM6K/sBUEQypIlS5bsUSF3ivabQp8xYwYaGxv7K3tBEISyhIi2hH0mLhdBEIQKQRS6IAiDhjW72vq7CCVFFLogCIOCpW/ux/m3PIOVO1r6uyglQxS6IAiDgl0tnQCAA4fK5a3U9ohCFwRhUNDakVPkqUy2n0tSOkShC4IwKGjtzCn0dKZyX3ciCl0QhEFBa0caAJDOioUuCIJQ1rgWekosdEEQhPLG9aGLhS4IglDmtHSIhS4IglARtHbmfOiZrCh0QRCEsibvcpFji4IgCOWNbIo6ENEFRLSGiNYT0fUhcc4momVEtJKIni5uMQVBEOIxGI4tFnx9LhElAdwK4DwA2wAsJqL7lVJveOKMBvAzABcopd4kogklKq8gCII13eksOlIZAGKhzwewXim1USnVDeBOAAsCca4AcLdS6k0AUEo1FbeYgiAI0Wnr7Hl/y2D/pugUAFs9/29zwrwcBWAMET1FREuI6EruQUR0FRE1ElFjc3NztBILgiBY4p5wASrb5WKi0IkJC05xVQBOAXAxgPMBfJ2IjtISKXW7UqpBKdVQX8/eoCQIglB03BMuQGW7XEyuoNsGYJrn/6kAdjBx9iilDgI4SETPADgBwNqilFIQBCEGrT6Xy+C20BcDmE1EM4moBsDlAO4PxLkPwNuIqIqIhgI4FcCq4hZVEAQhGu4JFwBIV/AXiwpa6EqpNBFdA+BRAEkAdyilVhLR1c7ntymlVhHRIwCWA8gC+JVSakUpCy4IgmCKz0KvYB+6icsFSqmHADwUCLst8P8PAPygeEUTBEEoDu3Opujw2qpBf8pFEAShrOlK586gD6tNVvSmqCh0QRAqnm5HiddVJyva5SIKXRCEiqc7nUVNVQJVCRKXiyAIQjmTymRRk0ygOpmQS6IFQRDKmbyFnqSKPrYoCl0QhIqnO52z0KsSYqELgiCUNalMFtVVhOqk+NCFAUzLoRTWN7X1dzEEYUDTlemx0OWUizBgue2ZDbjily/3dzEEYUCTSmdRncz50OUcujBgae1I+b7WLAiCTncmi9qq3CkXuSRaGLCkMtmK9gkKQjHwnkOXTVFhwJLOKKSzCkqJUheEMFKZnMulOpmQY4vCwKXbsTYq2S8oCHHxnUMXC10YqLjulkpeRgpCXLozCtXJBJIJ2RQVBjCuIhc/uiCE053OoKYqgWo5tigMZPIulwoWUkGIS3cmi9qk63KpXONHFHqZIy4XQShMKq3ym6KVPFaMFDoRXUBEa4hoPRFd30u8txBRhoguK14Rhd4Ql4sgFKY743l97mA+5UJESQC3ArgQwBwAHyaiOSHxvofc3aNCH5FyhLO7gq0OQYhLzzdFExVt/JhY6PMBrFdKbVRKdQO4E8ACJt5nANwFoKmI5RMKkEqLhS4IhehyLPTqJA36TdEpALZ6/t/mhOUhoikALgXguzg6CBFdRUSNRNTY3NxsW1aBIZU/h165QioIcVBKeb4pmkBWAdkKdbuYKHRiwoKtcQuA65RSmd4epJS6XSnVoJRqqK+vNyyi0BuuP1AUuiDwuGOkJkmoSubUWaWeCqsyiLMNwDTP/1MB7AjEaQBwJxEBwHgAFxFRWil1bzEKKYTT7bpcKtTiEIS4uGOkpqrHfk1nFGpNtF+ZYVKlxQBmE9FMANsBXA7gCm8EpdRM928i+i2AB0SZ9w2uP9D1pQuC4MddvVYnE3BfeVSpe04FFbpSKk1E1yB3eiUJ4A6l1Eoiutr5vFe/uVBa3K8xp8RCFwQWr4Xuvjp3MLtcoJR6CMBDgTBWkSul/i1+sQRTXMtcLHRB4OlyFXoygZSz/TdoLXRhYONaGpV8FEsQ4uC6XGqqElCo7FNhotDLnLzLpUItDkGIi/ulu5pkAlnHiV6phwhEoZcx2azq8QlWqMUhCHFJpXNjxHu5RaW+E10Uehnj3dipVJ+gIMSlO5P7ekxNVc+rcyvVQpe3LZYxXjeLvMtFEHi6PKdcqhI5lVepBpAo9DLGu2ys1CWkIMTFNXyqnfehA5V7bFEUehnjtcplU1QQeNxz6LVVufehA2KhCwMQrxKvVItDEOLi/aZoVSJnoVfqilYUehnjd7lUpsUhCHHxflO0yrHQK/Wb1aLQy5iUz+VSmRaHIMTFp9DFQhcGKj6Xi1jogsDSnXe5eF6fW6HjRRR6GSMWuiAUJr8pmkz2bIpW6J6TKPQyJiXHFgWhIPlN0SryuFzEQhcGGP5TLpUpoIIQl27P2xZ7LPTKHC+i0MsYn8tFXp8rCCzdmSwSBFR5vlhUqStaUehljHfZWKkWhyDEpTuTzVvm7lf/K3VFa6TQiegCIlpDROuJ6Hrm848Q0XLn5wUiOqH4RRWCeL8pKu9yEQSe7nQ2f5+o+7sr1et99mVLQYVOREkAtwK4EMAcAB8mojmBaJsAnKWUmgfgWwBuL3ZBBR3X5ZKgyl1CCkJcUpksahwLfWRdFWqqEmhu6+rnUpUGEwt9PoD1SqmNSqluAHcCWOCNoJR6QSm13/n3JQBTi1tMgcN1uQytqarYc7WCEBevhU5EmDSyDrtaO/u5VKXBRKFPAbDV8/82JyyMTwB4mPuAiK4iokYiamxubjYvpcDiulmG1CTlHLoghHCoO4O66mT+/0kj67CrZfAqdGLCWHOQiM5BTqFfx32ulLpdKdWglGqor683L6XA0mOhJyv2XK0gxGVnSycmjazL/z9xVB12D2ILfRuAaZ7/pwLYEYxERPMA/ArAAqXU3uIUT+gN1yofUi0WuiCEseNAByaPHpL/f9LIWuxq7YRSlWcEmSj0xQBmE9FMIqoBcDmA+70RiOhwAHcD+KhSam3xiylwuEp8aE2yYo9hCUIcUpksdrd2YsqYHoU+cWQdOlNZtHak+7FkpaHgnaJKqTQRXQPgUQBJAHcopVYS0dXO57cBuAHAOAA/IyIASCulGkpXbAHo+abo0Joq7DvY3c+lEYSBx66WTmQVMGV0j8tl0qjc37taOzFqaHV/Fa0kGF0SrZR6CMBDgbDbPH9/EsAni1s0oRApz6Zouk1cLoIQZPuBDgAIuFx6FPrRk0b0S7lKhXxTtIxJZ7IgAuqqk3JsURAYdjgKfcpov8sFAHZX4EkXUehlTHdGoTqZQHWSZFNUEBh2MBb6RI+FXmmIQi9jUpksqhOE6kRCji0KAsP2Ax0YP7zGdw69piqBccNq8sq+khCFXsakM1lUV+XeICcWuiDobD/Q6bPOXY6bMgqvvnmg7wtUYkShlzE9LpeEKHRBCKCUwoamdkwbM1T77LQjxmLN7jbsba+sd7qIQi9j0q7LJUmyKSoIAbbsPYTtBzpw2hFjtc9OO2IcAOCljfv6ulglRRR6GZNyXC7VyUTF3pEoCFF5dl3ufVFnztZfMzJ3yigMq0nipY2V9aV2UehlTCqjUJUgVCUTSGVURX6VWRCi8uy6PZg6ZghmjNNdLtXJBObPHIsnVu2uqFdPi0IvY1LOTSzV7sW38vV/QQAAdHRn8OKGvXjb7PFwvr2uccWp07GjpRMPvr6zj0tXOkShlzGpTO49z9XOu57l6KIg5Lhz8Zto60rj0pPCr2Y495gJOLJ+GG57eiMyFWIMiUIvY/IuF8dCl2voBAHoSmfwi6c3Yv7MsZg/U98QdUkkCJ89dzZW7WzFz59a34clLB2i0MuYQ91pVCcT+dtYKskXKAhR+c6Dq7CrtRPXnju7YNz3njAZ7z1hMn70+Do8uLz8XS+i0MuUQ91prNjeirlTRvXcZC4uF2GQ84cXN+N3L27BJ86cibfOGl8wPhHhO++bi5MPH43P/GUpfvrkurI2jEShlykvrN+L7kwW5xwzASPqci/N3H7gUD+XShD6h85UBt9+4A18/b6VeMcxE3DdBccYpx1eW4Xff/xUXDxvMm5+bC0u+PGzuHvpNhzqLr/3pRu9PlcYeCxa04RhNUk0zBiD7nQWQ2uS+OvirThlerjPUBAqjY7uDO5/bTt+umg9tu7rwEdPm44b3jMH1Uk7W3VITRI/+fBJuHjuYbj5sTX4z7+9hq/fuwLnHzcJb501HvNnjMW0sUNCT8wMFEShlyH7DnbjsTd244xZ41FblURtVRILTpyMe17djq9eNKfiXtovCC5KKWzb34HFm/fhidVNeGZNM9q60pg7ZRRu+uQ8nGHgZumNC46fhHfNmYhXNu/DPUu349E3duHuV7cDAEbUVuGICcMxq344Zo4fivoRtRg/3PkZUYsxQ6tRV5VEItF/Sp9MvoxCRBcA+DFyNxb9Sil1U+Bzcj6/CMAhAP+mlFra2zMbGhpUY2Nj1HIPWlbtbMUX//4a1je148//fhpOmT4GALByRwve/ZPncPLhY/Dzj5yMCZ5LcQWhnFBK4cChFLYf6MD2Ax3YcaAD2/Z3YNXOVqzc0YqWjhQAoH5ELc45uh7vO3kqTp05tiTWczarsK6pHa9s3od1u9uwobkdG5oO9vrq3SHVSQytSWJITe53TVUCyUQCSQKSCUKCCAtOnIIrTj08UpmIaEnYjXAFLXQiSgK4FcB5yF0YvZiI7ldKveGJdiGA2c7PqQB+7vwWIpDOZNHelcau1k7sasn9rGtqxyub9uH17S0YUVeF2/7llLwyB4DjJo/CrVecjGv/ugxnfm8Rzj66Hg0zxmDe1NGYPGoIJoys9b1CVBCiopRCKqOQySqkslmkMwrpTBaprEImo9CdyaIzlUFHKpP73Z37uyuVRYcT3tGdQVtnGgc6utHakcKBQykc6EihpSOFlkMp7QhuXXUCR00cgYvmTsKcyaNw4tTROG7yyJJbw4kE4ehJI7SbjTpTGexp78Ke9m7saetCc3sXWjpSONSdQUd32vmdwaHuDLrSGWRUbnLIZBUySkGhNAcYTFwu8wGsV0ptBAAiuhPAAgBehb4AwO9Vztx/iYhGE9FhSqminwNau7sNDy7fmWsO1dMsSiHfSLm/A+FOgALyX5Hn4nkXLPl4zue5v5WWDt58FZBROQFPZ1VO2LMKmWzP/+5AyOQ/zwn6wa4MDnal0ZXWd9lrqxKYO2UUrrvgGFwx/3DWrXLR3MMw57CR+M3zm/DU2mY89sZu3+ejhlRj1JBqDK1Jos61IqpzlkR1MoFkgpAkQjLp/E7kfqoShIT7mwiuIUTI/UEEuMMq/1nAWsrFIX+cXj7raVt/v6hePutJpwJx/PGDn7kB3LPd/8M+C5OrYFmD+WplBZBVClmVC8tmA/875cuqXHhYnJ7Pc3m6cZRy43rj+9O7spvJKqQy2bxsunKbcmS6WF/CGV5bhVFDqjF6aE4uj5o4HKOGVGPkkGpMGFGHKaOHYMroIZg8ug5jh9UMKP91XXUSU8cMxVTmTY79iYlCnwJgq+f/bdCtby7OFABFV+jrdrfjx0+sy//vVybEKgk4cbyKyKtUXEEhT9z880LiEeVTaAqtKpFTjlXJnBJMJnK3CrnKsSqRQG11Vf6zqgRhaE0SQ2uTGFZThaE1VRheV4WJI2tx2Kg6TBo1BBNH1KLKYKNnxvhh+OaC4wEAzW1dWLWzFbtbO9HU1oXdrZ1o7UihI5WzHDpTGbR0pNDRnclNMBmVH9C+H09YhXyhjiUoSwA/6QR+hcoV95ygfHnTJ4lAREgQkHB+ExESidyz3XDKf+75O+F+5qSDP04yQajOT8j+POB5diKBnLy68pt0/07k5bnKkeeqZE52q4OfVSVyhkJ1EnXVCdQ5RkOdEzakOonaqkS/+porFROFzrV6cFibxAERXQXgKgA4/PBo/qOL5k7C5psujpR2sFE/ohb1I/Q3zcWBt357rEz/Z7pV2rv12/McTiFyKwPk4/GfhU323mcLQqVgotC3AZjm+X8qgB0R4kApdTuA24HcpqhVSR1kEPYvQavSCe2XsgiC4MfksOZiALOJaCYR1QC4HMD9gTj3A7iScpwGoKUU/nNBEAQhnIIWulIqTUTXAHgUuWOLdyilVhLR1c7ntwF4CLkji+uRO7b4sdIVWRAEQeAwOodekoyJmgFsiZh8PIA9RSxOuSPt4UfaowdpCz+V0B7TlVLs5li/KfQ4EFFj2MH6wYi0hx9pjx6kLfxUenvIy7kEQRAqBFHogiAIFUK5KvTb+7sAAwxpDz/SHj1IW/ip6PYoSx+6IAiCoFOuFrogCIIQQBS6IAhChVB2Cp2ILiCiNUS0noiu7+/y9AVEdAcRNRHRCk/YWCJaSETrnN9jPJ99xWmfNUR0fv+UujQQ0TQiWkREq4hoJRF9zgkfdO1BRHVE9AoRvea0xTed8EHXFl6IKElErxLRA87/g6c9lPMqznL4Qe6bqhsAHAGgBsBrAOb0d7n6oN5vB3AygBWesO8DuN75+3oA33P+nuO0Sy2AmU57Jfu7DkVsi8MAnOz8PQLAWqfOg649kHuJznDn72oALwM4bTC2RaBd/hPAnwE84Pw/aNqj3Cz0/LvZlVLdANx3s1c0SqlnAOwLBC8A8Dvn798BuMQTfqdSqksptQm51zHM74ty9gVKqZ3KuQ1LKdUGYBVyr2oedO2hcrQ7/1Y7PwqDsC1ciGgqgIsB/MoTPGjao9wUeth71wcjE5XzAjTn9wQnfNC0ERHNAHAScpbpoGwPx72wDEATgIVKqUHbFg63APgyAO8tMYOmPcpNoRu9d32QMyjaiIiGA7gLwLVKqdbeojJhFdMeSqmMUupE5F5ZPZ+Iju8lekW3BRG9G0CTUmqJaRImrKzbo9wUutF71wcJu4noMABwfjc54RXfRkRUjZwy/5NS6m4neNC2BwAopQ4AeArABRi8bXEGgPcS0Wbk3LHvIKI/YhC1R7kpdJN3sw8W7gfwr87f/wrgPk/45URUS0Qzkbu4+5V+KF9JoNwNG78GsEop9UPPR4OuPYionohGO38PAfBOAKsxCNsCAJRSX1FKTVVKzUBONzyplPoXDKb26O9dWdsf5N67vha5Hemv9nd5+qjOf0HuftYUclbFJwCMA/AEgHXO77Ge+F912mcNgAv7u/xFboszkVsWLwewzPm5aDC2B4B5AF512mIFgBuc8EHXFkzbnI2eUy6Dpj3kq/+CIAgVQrm5XARBEIQQRKELgiBUCKLQBUEQKgRR6IIgCBWCKHRBEIQKQRS6UBEQUYaIlnl+ZhTpubcQ0dudv58iogbPZzO8b8Bk0tYT0SPFKIcgmFDV3wUQhCLRoXJfgddwvoxESqks93kYRDQWwGlKqWujFEgp1UxEO4noDKXU81GeIQg2iIUuVCSO9byKiH4GYCmAaUT0JSJaTETL3XeHO3G/6rwP+3Ei+gsRfdH56DIARhY2Ef3KszpoJqL/dj66F8BHilg1QQhFLHShUhjivHUQADYB+DyAowF8TCn1H0T0LuS+2j0fuZcy3e+4Ug4i9zXxk5AbD0sBuC93OgPAPwL5/ImIOpy/a+C81U8p9UkAIKLpAB4F8FsnTiOAbxetloLQC6LQhUrB53JxfOhblFIvOUHvcn5edf4fjpyCHwHgHqXUISed991AhwFoDuTzEaVUoyePBzx51gH4O4BrlFJbnOAmAJNj1k0QjBCFLlQyBz1/E4DvKqV+4Y1ARNci/JWpHQDqLPK7DcDdSqnHPWF1znMEoeSID10YLDwK4OPOe9RBRFOIaAKAZwBcSkRDiGgEgPd40qwCMMvk4UT0/wCMUErdFPjoKORenCUIJUcUujAoUEo9htw9ky8S0evI+cZHqNx1dn9F7q2NdwF41pPsQeTe2mfCFwHM9WyMXu2En+M8RxBKjrxtURA8ENE3ALQrpW52/n8OwLtV7gKJKM97BsACpdT+ohVSEEIQC10QeucLAA6PkpCI6gH8UJS50FeIhS4IglAhiIUuCIJQIYhCFwRBqBBEoQuCIFQIotAFQRAqBFHogiAIFcL/B7+3/8SvFVj0AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#导入使用的库\n",
    "import numpy as np\n",
    "import pylab as pl\n",
    "\n",
    "fs = 900#自己改频率\n",
    "N = 512\n",
    "t = np.arange(0, 1, 1/fs)\n",
    "x = np.sin(2*np.pi*100*t)+2*np.sin(2*np.pi*400*t)\n",
    "\n",
    "xs = x[:N]\n",
    "xf = np.fft.rfft(xs)/N\n",
    "freqs = np.linspace(0, int(fs/2), int(N/2+1))\n",
    "xfp = np.abs(xf)\n",
    "pl.figure()\n",
    "pl.subplot(211)\n",
    "pl.plot(t[:N], xs)\n",
    "pl.xlabel(\"Time(S)\")\n",
    "pl.subplot(212)\n",
    "pl.plot(freqs, xfp)\n",
    "pl.xlabel(\"Freq(Hz)\")\n",
    "pl.show()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "655bb755",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
